Capture of field stars by globular clusters in dense bulge regions 
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ABSTRACT 

The recent detection of a double Red Giant Branch in the optical color-magnitude 
diagram (CMD) of the bulge globular cluster HP1 (Ortolani et al. 1997), a more pop- 
ulated metal-poor steep one corresponding to the cluster itself, and another metal-rich 
curved, led us to explore in the present Letter the possibility of capture of field stars 
by a globular cluster orbiting in dense bulge regions over several gigayears. Analyti- 
cal arguments, as well as N-body calculationsLatex for a cluster model of 10 5 M Q in a 
bulge-like environment, suggest that a significant fraction of cluster stars may consist of 
captures. Metal-poor globular clusters in the inner bulge, like HP1, contrasting at least 
in A[Fe/H] = 1.0 dex with respect to the surrounding metal-rich stars, are ideal probes 
to further test the capture scenario. In turn, if this scenario is confirmed, the double 
RGB of HP1 could provide direct estimates of blanketing amounts, which is fundamental 
for the photometric calibration of metal-rich stellar populations. 

Subject headings: Galaxy: kinematics and dynamics — Galaxy: globular clusters: indi- 
vidual (HP1) — Galaxy: globular clusters: general — methods: analytical — methods: 
numerical 
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1. Introduction 

Recently, Ortolani, Bica & Barbuy (1997, hereafter 
OBB97) studied V vs V-I diagrams of the inner bulge 
globular cluster HP 1 , which revealed the presence of 
a double Red Giant Branch (RGB). One RGB is more 
populated and steep, and is accompanied by a well- 
developed Blue Horizontal Branch (HB), which un- 
doubtedly characterizes HP 1 as a metal-poor globu- 
lar cluster. The other RGB, observed in a r< 23" 
extraction, is curved and extended, characteristic of 
a nearly-solar metallicity stellar population, like that 
of the globular cluster NGC 6553 or the Baade Win- 
dow itself (Ortolani et al. 1995). OBB97 favored the 
interpretation of a bulge contamination over that of a 
merger of two clusters, to explain the secondary RGB 
in HP 1. In the present Letter we explore a scenario 
where this contamination may correspond to physical 
captures of bulge stars. 

Van den Bergh (1996) proposed the merger sce- 
nario to explain composite (blue and red) Horizontal 
Branches for some intermediate and low metallicity 
globular clusters in the Milky Way. For such mergers 
to occur, individual clusters must have low relative 
velocities. This condition is matched in dwarf galax- 
ies, where the clusters could have merged in a first 
step inside the dwarf, and subsequently would have 
been accreted by our Galaxy. As pointed out by van 
den Bergh (1996) the exceptionally bright globular 
cluster M 54 in the Sagittarius dwarf (Ibata, Gilmore 
& Irwin 1994) might be such an example. 

The merger scenario would not be appropriate for 
the particular case of HP 1 since (i) the bulge veloc- 
ity dispersion is very high, which makes unprobable 
mergers of globular clusters within the bulge, and (ii) 
the stellar population of the secondary RGB in HP 1 
is exceedingly metal rich to have been originated in a 
dwarf galaxy. 

In the present Letter we study, both analytically 
and with N-body simulations, the possibility of cap- 
ture of bulge stars by a globular cluster, as an al- 
ternative mechanism to explain the secondary RGB 
observed in the color-magnitude diagram (CMD) of 
HP1. 

2. Analytical approach 

We study the capture scenario by a globular clus- 
ter orbiting inside the bulge, assuming for the bulge 
total mass and mass distribution those predicted in 



the Galaxy model by Hcrnquist(1993; see also, Hern- 
quist 1990a). OBB97 derived for HP 1 a distance of 
~1.3 kpc from the Galactic center. For this distance 
the cluster faces a star density of O.33M0/pc 3 (Hern- 
quist 1993). For the remaining cluster parameters we 
consider a grid of values which encompasses probable 
properties for a bulge globular cluster. The rotation 
of the galactic bulge (Menzies 1990) implies that glob- 
ular clusters in the bulge might present considerable 
amounts of streaming motion. So we consider three 
possible values for the cluster velocity (V c ): 50, 100 
and 160 km s -1 . For cluster diameters we test 10, 15 
and 20 pc, whereas for masses the values are 10 4 , 10 5 
and 10 6 M Q . For these ranges of sizes and velocities 
the cluster encounters a total number of bulge stars 
in the range 1.3 x 10 3 < N b < 7.1 x 10 4 during 1 
Myr (assuming 1 M for bulge stars). Similarly, the 
size and mass ranges imply cluster escape velocities 
3kms _1 < V e < 42kms~ 1 . The cluster mass is the 
most important parameter determining the escape ve- 
locity V e , for the adopted grid. 

The bulge velocity distribution is assumed to be 
isothermal (Binney & Tremaine 1987) with a disper- 
sion of mkms- 1 (e.g. Sharpies et al. 1991). The 
number of bulge stars capable to be captured by 
the cluster must satisfy the following velocity con- 
straints: \V C \ - \V e \ < \V S \ x cos < \V C \ + \V e \ and 
\V S \ x sin 9 < \V e \ where 6 is the angle between the 
star velocity (V s ) and that of the cluster (V c ). Adopt- 
ing a typical mass of 1 M Q for bulge stars, the number 
of stars N s which are captured over 1 Gyr by the clus- 
ter can be summarized as follows: (i) for a 1O 6 M 
cluster with escape velocity V e = 34kms _1 and 
streaming velocity V c = 50kms _1 , N s = 1.2 x 10 5 , 
whereas for V c = 100 kms" 1 , N s = 4.1 x 10 4 . 
(ii) For a 10 5 M Q cluster with escape velocity V e = 
llkms -1 and streaming velocity V c = 50kms _1 , 
N s = 5.1 x 10 3 , whereas for V c = 100 km s" 1 the 
captures become negligible. Finally, a 10 4 M Q cluster 
docs not capture a significant number of stars. This 
would set constraints on the mass of a cluster with a 
secondary RGB, like HP 1. 

3. Numerical simulation 

In order to further check the capture scenario, 
we performed a N-body numerical experiment. We 
used a hierarchical tree algorithm (Barnes & Hut 
1986) in the CESUP-UFRGS Cray Y-MP2E com- 
puter. The algorithm was optimized for vector archi- 
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Fig. 1. — Initial distribution of cluster and bulge particles encountered during 10 Myr, respectively separated in (a) 
and (b) for clarity purposes, (c) and (d) show respectively the cluster and nine bulge particles (90 M Q ) captured 
after 10 Myr. Finally, (e) and (f) test the cluster and capture stability for further 40 Myr. 
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tectures (Hernquist 1987,1990b). The adopted toler- 
ance parameter is 0.7, and the calculation includes 
quadrupole terms. 

The globular cluster was modeled with a Plum- 
mer polytrope with 10 4 particles amounting a mass of 
10 5 M Q . The Plummer cluster cutoff diameter is 15 pc 
and the core diameter is 2pc. The cluster was left 
to evolve isolated during 1 relaxation time (2.8 Myr) 
to check its stability. Due to computational limita- 
tions on the number of particles in the simulation, we 
adopted for cluster and field stars 10 M Q particles. 
More details on similar star cluster simulations are 
given in Rodrigues et a/.(1994). 

The adopted Galactic bulge density at the cluster 
location is that of the previous section. A cluster with 
a streaming velocity of 50kms _1 encounters a mass 
of bulge stars of 2948 M Myr _1 . The bulge stars 
encountered by the cluster during 10 Myr were placed 
in a homogeneous sphere with size equal to that of the 
cluster, as shown in Figs. [l]a,b. The speed |V^| of a 
given star in the bulge sphere was randomly sampled 
from a Gaussian distribution with a = 113 km s -1 
(Sect. ||). The Cartesian velocity components were 
calculated from \V S \ assuming isotropy. 

The simulation begins with spatially coincident 
cluster and bulge spheres, placed to orbit around a 
massive particle of (7 x 1O 9 M ), representing the 
bulge mass internal to the cluster position (Sect. ||). 
The status of the simulation after 10 Myr is shown in 
Figs. 0c, d, respectively for the cluster and nine bulge 
particles (equivalent to 90 M Q ) which were captured. 
In order to further check the capture stability we left 
the simulation evolve up to 50 Myr (Figs. |l|c,d) and 
seven particles (70 M ) remained trapped. We re- 
mind that during these supplementary 40 Myr no new 
encounters with bulge stars were considered. The cap- 
ture rate obtained from this numerical simulation im- 
plies that about 7000 to 9000 M can be captured by 
a 10 5 M Q globular cluster lurking in the bulge during 
1 Gyr, in good agreement with the analytical results 
of Sect. |. 

4. Discussion and concluding remarks 

In the present Letter we proposed a capture mech- 
anism to explain a metal-rich secondary RGB in the 
CMD of the bulge metal-poor globular cluster HP 1. 
So far the proposed mechanisms affecting the number 
of stars in globular clusters are the escape of stars by 
evaporation (e.g. Spitzer & Thuan 1972) and ejection 



(e.g Henon 1969). By far evaporation is the most im- 
portant escape mechanism: typically, a globular clus- 
ter with 10 6 Mq should be stable during a Hubble 
time, whereas one with 10 5 M Q should lose an im- 
portant fraction of its mass. The disk shocking pro- 
cess (Ostriker et al.1972; see also Binney & Tremaine 
1987), not only by the disk but also by the bulge itself 
in such central regions, is another possible loss mech- 
anism. The presently proposed accretion mechanism 
for bulge globular clusters would balance such losses, 
allowing those with ~ 10 5 M Q to survive in the bulge. 
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Fig. 2. — Color-magnitude diagram of HP 1 from 
OBB97, where in addition to the mean locus of 
NGC6752 (dashed line), we superimposed that of 
NGC 6553 (continuous line) which fits the metal-rich 
sequences 

Depending on the cluster mass, as well as on the 
evaporation and capture rates it is possible to envis- 
age changes of the cluster stellar population content 
over a Hubble time. A very massive cluster, initially 
metal-poor, would conserve this character and add a 
secondary metal- rich component to its CMD, as could 
be the case of HP 1 . On the other hand, a less mas- 
sive initially metal-poor cluster might recycle its stel- 
lar content, if capture dominates over evaporation ef- 
fects, so that its CMD would become similar to that 
of the bulge. In this case the metal-poor evolution- 
ary sequences like blue HB and vertical RGB would 
become relatively less populated. Hubble Space Tele- 
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scope color-magnitude diagrams of metal-rich bulge 
clusters would be important to check this scenario. 

The metallicity distribution of the galactic globu- 
lar clusters is skewed towards high metallicities, or it 
is perhaps bimodal (Zinn 1980), and cannot be de- 
scribed by a one-zone metal enrichment model (Bica 
& Pastoriza 1983). The above described scenario of 
bulge star captures by globular clusters and recycling 
of their stellar population provides a natural way to 
create an asymmetrical histogram of cluster metal- 
licities from an initially Gaussian distribution. It is 
possible to speculate the existence of two families of 
metallic clusters: (i) genuine ones formed from en- 
riched gas, and (ii) recycled ones which were initially 
metal-poor and, over a Hubble time, developed CMD 
sequences like those of the Baade window by means 
of captures. 

The relative loci of evolutionary sequences in the 
HR diagram for globular clusters of different metallic- 
ities is fundamental for the calibration of stellar pop- 
ulation parameters. Da Costa & Armandroff (1990) 
studied metal-poor globular clusters together with 
47Tucanae, whereas Bica, Barbuy & Ortolani (1991) 
studied those of nearly solar metallicity. The HP 1 
CMD containing combined metal-poor and metal-rich 
stellar populations could provide a means to directly 
measure relative blanketings, independently of stel- 
lar atmosphere and spectral models. We show in 
Fig. | the CMD of the spatial extraction for r<23" 
from OBB97, where in addition to the mean lo- 
cus of a metal-poor globular cluster (NGC 6752 with 
[Z/Z Q ] =-1.54) fitted to the sequences of HP 1 itself, 
we also superimpose that of the nearly-solar metallic- 
ity globular cluster NGC 6553 (Ortolani et al. 1995). 
This HP1 CMD was built from images obtained at 
the ESO NTT 3.55 m telescope and was reduced with 
the DAOPHOT II package with particular care for 
crowded field extractions and calibrations (Ortolani 
et qZ.1996, OBB97). 

We illustrate relative measures of reference points 
in the CMD. The magnitude difference between the 
metal-poor and metal- rich horizontal branches is AVhb ~ 
0.55 mag, which has important implications for the 
relative distances of globular clusters (e.g. Jones 
et aZ.1992). The magnitude difference between the 
brightest giants in the metal-poor and metal-rich 
RGBs is AVbg ~ 1.75 mag. However, if the tips of 
the metal-poor and metal-rich giant branches corre- 
spond to a similar stellar temperature, the difference 
could be as large as AVtip ~ 3.6 mag, due to blan- 



keting effects. Finally, the color difference could be as 
large as A(V — I)tip ~ 2.2 mag. Such magnitude and 
color differences are fundamental quantities for blan- 
keting calibrations of metal-rich stellar populations. 
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